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Abstract 

We study the existence of particular traveling wave solutions of a nonlinear parabolic 
degenerate diffusion equation with a shear flow. Under some assumptions we prove 
that such solutions exist at least for propagation speeds c e]c*,+oo[, where c* > is 
explicitly computed but may not be optimal. We also prove that a free boundary hy- 
persurface separates a region where u = and a region where u > 0, and that this free 
boundary can be globally parametrized as a Lipschitz continuous graph under some 
additional non-degeneracy hypothesis; we investigate solutions which are, in the region 
u > 0, planar and linear at infinity in the propagation direction, with slope equal to 
the propagation speed. 



1 Introduction 

Consider the advection-diffusion equation 

d t T -V ■ (XVT) + V ■ (VT) = 0, (t,X) eR + xR d (1.1) 

where T > is temperature, A > is a diffusion coefficient and V = V(x\, ...,Xd) G K d is a 
prescribed flow. In the context of high temperature hydrodynamics, the diffusion coefficient 
A cannot be assumed to be constant as for the usual heat equation, but rather of the form 
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A = A(T) = X T m for some conductivity exponent m > depending on the model, see \T7\. 
We will consider here the case m/1. In Physics of Plasmas and particularly in the context 
of Inertial Confinement Fusion, the dominant mechanism of heat transfer is the so-called 
electronic Spitzer heat diffusivity, corresponding to m = 5/2 in the formula above. 

Suitably rescaling one may set A = m + 1, yielding the nonlinear parabolic equation 

d t T — A (T m+1 ) + V ■ (VT) = 0. (1.2) 

When temperature takes negligible values, say T = e — > 0, then the diffusion coeffi- 
cient A(T) = X T m may vanish and the equation becomes degenerate. As a result free 
boundaries may arise. We are interested here in traveling waves with such free boundaries 
T = d{T > 0} 7^ 0, and in addition T — » +oo in the propagation direction. 



When V = (1.2) is usually called the Porous Medium Equation - PME in short - 

d t T - A (T m+1 ) = (PME) 

and has been widely studied in the literature. We refer the reader to the book [16] for 
general references on this topic and to [U [21 E] for well-posedness of the Cauchy problem and 
regularity questions. As for most of the free boundary scenarios, we do not expect smooth 
solutions to exist, since along the free boundary a gradient discontinuity may occur: a main 
difficulty is to develop a suitable notion of viscosity and/or weak solutions; see [13] for a 
general theory of viscosity solutions and [6] in the particular case of the PME, [16] for weak 
solutions. 

The question of parametrization, time evolution and regularity of the free boundary for 



(PME) is not trivial. It has been studied in detail in [HI El [10] • When the flow is potential 



V = V$ (1.2) has recently been studied in [15], where the authors investigate the long time 



asymptotics of the free boundary for compactly supported solutions. 

We consider here a two-dimensional periodic incompressible shear flow 

V(x,y)= ^, a(y + l) = a(y) 

i 

for a sufficiently smooth a(y), which we normalize to be mean-zero J a(y)dy = 0. In this 

o 



setting (1.2) becomes the the following advection-diffusion equation 

d t T - A(T m+1 ) + a{y)d x T = (AD-E) 



with 1-periodic boundary conditions in the y direction. 



For physically relevant temperature T > it is standard to use the so-called pressure 
variable u = 2^Lj' m which satisfies 

m ' 

dtu — muAu + a(y)d x u = \ Vm| 2 . (1-3) 
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Remark 1. When m = 1 pressure u is proportional to temperature T , and this particular 
case will not be considered here. 



Looking for wave solutions u(t, x, y) = p(x + ct, y) yields the following equation for the 
wave profiles p(x, y) 

-mpAp+(c + a)p x = \\7p\ 2 , (x,y)eRxT l . (1.4) 



In the case of a trivial flow a = it is well-known [16] that for any prescribed propagation 
speed c > there exists a corresponding planar viscosity solution given by 

p(x,y) = p c (x) = c[x] + , (1.5) 

where [.]+ denotes the positive part. This profile is trivial for x < and linear for x > 0, 
with slope exactly equal to the speed c. The free boundary T = {x = 0} moves in the 
original frame with constant speed x{t) = —ct + est, and the slope at infinity therefore fully 
determines the propagation. 

In this particular case the free boundary is non-degenerate, Vp = (c, 0) ^ in the "hot" 
region p > 0. The differential equation satisfied by the free boundary T in the general 
case was specified in [TO], where the authors also show that if the initial free boundary is 
non-degenerate then it starts to move immediately with normal velocity V = — Vp| r . 



In presence of a nontrivial flow a ^ a natural question to ask is whether (AD-E) 



can be considered as a perturbation of (PME). More specifically we are interested here in 



the following questions: (i) do y-periodic traveling waves behaving linearly at infinity and 
possessing free boundaries still exist? (ii) If so for which propagation speeds c > 0, and is 
it still true that the slope at infinity equals the speed c? (iii) Is the free boundary wrinkled 
by the flow and how can we parametrize it? (iv) Is the free boundary non-degenerate and 
what is its regularity? 

We answer the first three questions above, at least for propagation speeds large enough. 
Main Theorem 1. Let c* := — mina > 0: for any c > c* there exists a nontrivial traveling 



wave profile, which is a continuous viscosity solution p(x,y) > of (1.4) on the infinite 
cylinder. This profile satisfies 

1. if D + := {p > 0} denotes the positive set, we have D + ^ and p\ D + G C°°(D + ), 

2. p is globally Lipschitz, 

3. p is planar and linear p(x,y) ~ cx uniformly in y when x — > +oo. 

Moreover there exists a free boundary V = d(D + ) ^ which can be parametrized as follows: 
there exists an upper semi- continuous function I(y) such that p(x,y) > <^ x > I(y). 
Further: 

• If y is a continuity point of I then T fl {y = y } = (I(yo),yo)- 

• If yo is a discontinuity point then T fl {y = yo} = [L(yo), I(yo)] x {y — Ho}, where 

I(y ) := liminfJ(y). 

y->vo 
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The question 4 is sill open. The non- degeneracy of the pressure at the free boundary 
Vp| r 7^ and the free boundary regularity are closely related. For the PME it was dis- 
cussed in [HI EJ [in] • We have, however, a partial answer under some strong non-degeneracy 
assumption. 

Proposition 1.1. With the same hypotheses as in Theorem^ assume that the non- degeneracy 
condition p x > a > holds in D + for some constant a. Then the function I(y) defined zn[7] 

is Lipschitz, and T = d{p > 0} = y), x = 

Remark 2. The condition of linear growth at infinity is natural because it mimics the planar 



traveling wave (1.5) for the PME. Let us also point out, a posteriori, that this linearity 



appears very naturally in our proof , see Section^ 

Recalling that we normalized J a(y)dy = ,we will always assume in the following that 

T 1 

the propagation speed c > is large enough such that 

0<co<c + a<d (1.6) 

for some constants cq,ci. This is indeed consistent with c > c* = — mina > in the main 
Theorem [U 

The method of proof of Theorem [I] is standard. We refer the reader to [4| for a general 
review of this method and to [H] for the special case of the PME. The proof relies on a 
simple observation: if p > 5 > ( |1.4[ ) is uniformly elliptic; we shall refer in the sequel to any 
solution p > 5 > as a 5-solution. The main steps are the following. 



We first regularize (1.4) by considering its (^-solutions with (5 C 1 on finite cylinders 
\—L,L\ x T 1 , L ^> 1 with suitable boundary conditions. In Section [2] we solve this regular- 
ized uniformly elliptic problem, and derive monotonicity estimates of its solution in the x 
direction. In Section [3] we take the limit L — > +00 for fixed 5 > 0, and establish 

Theorem 1.1. For any 5 > small enough there exists a smooth 5 -solution p > 5 on the 
infinite cylinder such that lim p(x,y) = 5 and p(x,y) ~ cx uniformly in y. 

The linear behavior will be actually proved in section [5] for the final viscosity solution, but 
the proof can be easily adapted for the 5-solutions. We complete the proof of parts 1. and 2. 
of Theorem [l] in Section [4] by taking the degenerate limit 5 — > + . Section [4] also contains the 
analysis of the free boundary and the proof of Proposition LI In Section [5] we investigate 



the linear growth and planar behavior at infinity. We refine part 3. of Theorem [T] as 

Theorem 1.2. Both for the viscosity solution of Main Theorem [I] and the 5-solution of 
Theorem \l.l\ the following holds when x — > +00: 

1. p(x,y) ~ cx, p x (x,y) ~ c and p y (x,y) — > uniformly in y 

2. If 1 < m ^ N and N := [m], there exist q\...qN and q* G M. such that 

p(x, y) = cx + x (q\X~™ + ... + qN%~™^ + q* + o(l) 



4 



The second part is novel compared to the PME, for which the planar wave is exactly linear 
p = cx at infinity. Once again, we prove this statement for the final viscosity solution, but 
the proof extends to the (^-solutions. Section [6] is finally devoted to uniqueness of the wave 
profiles, and we establish 

Theorem 1.3. The 5-solutions of Theorem \l . 1\ are unique up to x -translations. 



2 ^-solutions on finite domains 



Here we solve (1.4) on truncated cylinders Dl = [—L,L] x T 1 , L ^> 1 with a uniform 



ellipticity condition p > 5 > 0. We show in this section that this ellipticity can be achieved 
by setting suitable boundary conditions, and we therefore consider the following problem 



< 5 < A < R 



—mpAp + (c + ct)p x 
p = A, 
p = B, 



|Vpp 



(x = — L), 
{x = +L), 



(2.1) 



where the constants A and B are specified later. 



We will show that any solution of (2.1) must satisfy p x > 0, and therefore p > A > on 
Dl. Thus (2.1) is uniformly elliptic. We prove this x-monotonicity of p using the following 



non-linear comparison principle, which relies on the celebrated Sliding Method of Berestycki 
and Nirenberg |oJ. 

Let a < b, Q =]a,6[xT 1 and for any function / e C 2 (Q) fl C(O) define the nonlinear 
differential operator 

$(/):= -m/A/ + (c + a)f x - | V/| 2 . (2.2) 
Theorem 2.1. (Comparison Principle) Ifu,v G C 2 (Q) nC(fi) satisfy u, v > in Q and 

u(a,y) < u(x,y) < u(b,y) 



V(x,y) en 



v(a,y) < v(x,y) < v(b,y) 



(2.3) 



then 



> $(u) (fi) 
M > v (dQ) 



mm m > max v 

x=b x=a 



U>V (H) . 



The proof is easily adapted from [5]: the fact that the equation is invariant under x- 
translations and that the domain is convex in this direction allows to compare translates of 
u and v. 



Condition (2.3) may seem quite restrictive at first glance, as it requires u, v to lie strictly 



between their boundary values: the following proposition ensures that this holds for any 



positive classical solution of problem (2.1). 



Proposition 2.1. Any positive solution p G C 2 (Dl) fl C(Dl) of problem (2.1) satisfies 

\/(x,y)eD L p(-L,y) <p(x,y) <p(L,y). 
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Proof. Assume that p is such a solution: since p > on the (compact) cylinder [-L, L] xT 1 , 
equation —mpAp + (c + a)p x — \ Vp\ 2 = can be considered as a uniformly elliptic equation 
Lp = with no zero-th order term: the classical weak Maximum Principle therefore implies 
on Dr that A = min p < p < max p = B, and the classical strong Maximum Principle 

ensures that the inequalities are strict m. D^. □ 
Corollary 2.1. There exists at most one positive solution p G C 2 (Dl) fl C(-Dl) of problem 



(2.1) 



Proof. Assume p\ ^ pi are two different solutions: $(pi) = &{jp2) = 0, min pi = B > A 



max pj and by previous proposition pi,p2 satisfy condition (2.3): Theorem 2.1 yields pi > pj 

x=a ' " 

and therefore p\= P2- □ 

We recall that a function p + G C 2 {VL) PlC(O) (resp. p~) is a supersolution (resp. subsolu- 
tion) if $(p + ) > (resp. $(p~) < 0), and construct below two different types of planar sub 
and supersolutions. 

An elementary computation shows that a planar affine function p + (x,y) = A + x + B + 
is a supersolution (resp. p~(x,y) = A~x + B~ is a subsolution) if and only if + (c + 



a)A + > (A + ) 2 (resp. A~,<). Due to hypothesis (1.6) this condition is satisfied as soon 
as < A + < Co (resp. A~ > c\ or A~ < 0): any affine function with positive slope A + < cq 
(resp. A" > ci) is hence a supersolution (resp. subsolution). 

We will also use some additional planar sub and supersolutions. For any xq G R, M > 
5 > and C > consider the following boundary value problem 



-muu" + Cm' = (u') 2 
Uc{x) : ^ m(— oo) = 

u{x ) = M. 



j i 



Using ODE techniques one easily shows that there exists a unique solution uc, which satisfies 
uc > 5, C > u' c > and u" c > for x G R. Defining p + (a;, y) := mc(i) one easily computes 
for < C < c (hence c + a(y) > C) 

$(p+) > - mMcU £ + - ( M ' c ) 2 = 0, 

and p + is therefore a planar supersolution. The same computation shows that if c + a(y) < 
C\<C then p~(x,y) = uc{x) is a planar subsolution $(j5 _ ) < 0. 



This allows us to build planar sub and supersolutions tailored to (2.1) as follows. Let 
5 > be a small elliptic regularization parameter, and define 

{—muu" + Cqu' = (u') 2 , 

u(-oo) = S, (2.4) 

u(0) = 1. 

If 

B:=p + (L), (2.5) 
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/\2 



similarly define 

{—muu" + C]V = [W , 
u(-oo) = 5, (2.6) 
u(L) = B. 

As discussed above p~ ,p + are a planar sub and supersolution on Dl = [-L, +L] x T 1 . They 
satisfy all the hypotheses of our Comparison Theorem 2.1 and therefore pr < p + . 
We prove below that, choosing 



A :-- 



P 



-L)+p-(-L) 



(2.7) 



there exists at least one solution p of problem (2.1) lying between p and p 



Theorem 2.2. Fix 5 > small enough: for L > /arge enough and A, B defined by 
(2.7)-(2.5) there exists a unique positive classical solution p £ C 2 (D L ) fl C 1 (D L ) of (2.1). 
Moreover, it satisfies p~{x) < p(x,y) < p + (x) on Dl and p £ C°° (-Dl)- 



Proof. Uniqueness is given by corollary 2.1 It was shown in [TT] that if there exist strict sub 
and sup ersolut ions p~ < p + then there is a classical solution p satisfying p~ < p < p + . Note, 
however, that we set C = Cq,Ci in (2.4)-(2.6). This corresponds to non-strict inequalities 
c &(p + ) > and $(p _ ) < 0, meaning that these particular sub and super solutions are not 
strict ones. It is however easy to approximate p by strict sub and supersolutions pf, where 
pf > p + and p~ < p~ are such that pf — > p uniformly on Dl when e — > + ; this can be 



done setting cq — e instead of Co and ci + e instead of c\ in (2.4)-(2.6) (and also approximating 
boundary conditions). 

All the hypotheses of Theorem 1 in [TT] are here easily checked, and we conclude that 
there exists at least one solution p e £ C 2,a (D L ) fl C l (D L ) such that p~ <p £ < pf on D L and 



satisfying boundary conditions p £ (—L, y) = A, p £ (+L, y) = B. By uniqueness (corollary |2. 1 ) 
this solution is independent of e, i-e p £ = p; taking the limit e — > + yields p~ < p < p + on 
Dl, and p is smooth on by standard elliptic regularity. □ 

As we let L — > oo in the next section, we need monotonicity of p in the x direction, as 
well as an estimate on p x uniformly in L, the size of the cylinder Dl- 



Proposition 2.2. The solution p(x,y) of (2.1) satisfies 



<Px<Ci 



(2.8) 



on Dl, where c\ > is given in (1.6) 



Proof, p £ C°°(Dl) DC 1 (D L ) is smooth enough to differentiate (1.4) with respect to x, and 
q=:p x E C°°(Dl) n C (D~l~) satisfies 



mpAq + [(c + a)q x — 2Vp ■ Vg] — (mAp)q = 0. 



(2.9) 



We first prove the upper estimate p x < c\. Since we set p(—L,y) = A < B = p(L,y) 
there exists at least a point inside Dl where p x > 0; any potential maximum interior point 
for q = p x therefore satisfies q > 0, and of course Vg = 0, Ag < 0. Using (2.9) we compute 
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at such a positive interior maximum point (mAp)q = —mpAq > 0, hence mAp > and 



mpAp > 0. The original equation (1.4) satisfied by p > yields now 

< mpAp = (c + a)p x - |Vp| 2 =^> (p x ) 2 < \ Vp\ 2 < (c + a)p x , 

and since q = p x > at this maximum point we obtain q = p x < (c + a) < c\. 

We just controlled any potential maximum value for p x inside the cylinder, and we control 
next p x on the boundaries using sub and supersolutions as barriers for p. Recall that the 
boundary values are flat, p(—L, y) = A and p(L, y) = B. 

On the right side we use the previous subsolution p~ (x) as a barrier from below: recalling 
that p~ < ci and that p~ and p agree at x = L we obtain p x (L, y) < p x (L) < c\. 

On the left we use a new planar supersolution as a barrier from above: let p(x) be the 
unique afline function connecting p(—L) = A and p(L) = B, with slope s = ^wA- Using 



(2.4)-(2.6) it is easy to see that B ~ cqL and A ~ 5 when L — > +oo for fixed 5. As a 



consequence s ~ cq/2 < cq for L large enough, and p is indeed a supersolution. Since p 



agrees with p on both boundaries our comparison Theorem 2A_ ensures that p < p on Dl, 
thus p x < s < c < Ci on the left boundary. 

The lower estimate q = p x > is inside D L & classical consequence of the Sliding 
Method [S]. In order to establish this strict monotonicity up to the boundaries we consider 
—mpAp+ (c + a)p x — \ Vp\ 2 = as a linear elliptic equation Lp = with trivial zero-th order 



coefficient. Proposition |2. 1| and flatness of the boundaries show in addition that p attains a 
strict minimum at any point of the left boundary, and also a strict maximum at any point 
of the right boundary. Hopf Lemma then implies that p x > on both boundaries. 

□ 

Proposition 2.3. (Uniform Pinning) There exists large constants K > 0, Ki m K — \/K 

and K 2 ~ K + \/K such that, for any L large and any 5 small enough, there exists x* = 
x*(L,5) G]0, L[ such that 

1. lim (L — x*) = +oo, 

L— >+oo 

2. K x < v (x*,y) < K 2 . 



The constants K, Ki,K 2 depend on the upper bound c\ in (1.6) and m > 0, but not on L or 
5. 

Proof. The idea is as follows. When x increases from — L to L the map x h > J p(x, y)dy 

T 1 

increases from y4~5<ltoi?~ c L ^> 1. For fixed large K and L large enough this integral 
has to take on the value K at least for some x g] — L,L[. The equation for p then allows 
us to control the ?/-oscillations of p along this line by 0(\/K). If K is chosen large enough 
these oscillations are small compared to the average, and p along this line will therefore be 
of the same order O(K) than its average J pdy. This will be our pinning line x = x* (up to 
some further small translation). 

• Choose a large constant K > 1, and for x G [— L, L] define F(x) := f p(x,y)dy: since 

p(0,y) < P + (0) = 1 we have that -F(O) < K. By convexity p~ lies above its tangent 
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plane ti{x) at x = L, and we recall that we had set p (L) = p(L, y) = p + (L) = B. For 
L large and 5 small tAx) = K has a unique solution x = xk given by xk = L + K ~,f s , 
and F(xk) > P~0ekt) > ^l{xk) = -K"- Remarking that F is increasing {p x > 0), there 
exists a unique x^(L,5) g]0,x^] such that 

F(x* K ) = J p(x* K ,y)dy = K. 



Manipulating (2.4)-(2.6) it is easy to check that for K,5 fixed and L — > +00 there 
holds 

B = p + (L) ~ c L 1 T K-B / c 

PxW ~ ci J V ci 

as a consequence of (1.6) the line (5, L) stays away from both boundaries. 

Let us now slide the whole picture to the left by setting p(x,y) = p(x + x* K ,y), so 
that corresponds in this new frame to x = 0; for simplicity of notation we 

will use p(x, y) instead of p(x, y) below. The corresponding domain still grows in both 
directions when L — > +00, and by definition of x* K we have that 

p(Q,y)dy = K. 



T 1 



We claim now that there exists a constant C, depending only onm^l and the upper 
bound for the flow ci, such that 

Vx>0, J J \Vp\ 2 dxdy < C(K + x). (2.10) 

[0,1] xT 1 

Indeed, integrating by parts the Laplacian term in —mpAp + (c + a)p x = | Vp| 2 over 
a subdomain Q = [0,x] x T 1 and combining the resulting |Vp| 2 term with the one on 
the right hand side yields 

(m — 1) J J I Vp\ 2 dxdy + m J pp x (0, y)dy — m J pp x (x, y)dy + J J (c + a)p x dxdy = 0. 
n t 1 t 1 n 

(2.11) 



1. If m — 1 > we use m J pp x (0, y)dy > and J J (c + a)p x dxdy > in (2.11 ). This 

t 1 n 
leads to (m — 1) ff \Vp\ 2 dxdy <mf pp x (x,y)dy, and since < p x < c\ 



/ / \Vp\ 2 dxdy < — — - / p(x,y)dy. 



T 1 



Our monotonicity estimate < p x < c\ again yields 

p(x, y)dy = J p(0, y)dy + J J p x dxdy < K + a 



x. 



and together with the previous inequality 

Vx>0, // \Vp\ 2 dxdy < -^-(K + Cl x) < C(K + x) 
J J m — 1 

[O^xT 1 



2. If < m < 1 we use pp x (x,y) > in (2.11) to obtain 



(1 — m) J J |Vp| dxdy < m J pp x (0, y)dy + J J (c + a)p x dxdy. 
n t 1 n 

Since < p x < C\ and < c + a < C\ this leads to 

J I \Vp\ 2 dxdy < f^J p(0, y)dy + J I c\dxdy 
n t 1 n 

with C = max(mci, c^) depending only on m and C\. 



max y) — min p(x, y) 
in the y direction: by Cauchy-Schwarz inequality we have that 

2 



In the spirit of [12] we control now the oscillations 0(x) 
: t 
/ 



2 (x) < 



VET 



< / \Py(x,y)\ 2 dy < / |Vp| 2 (x, y)dy 



T 1 



T 1 



and integrating from x = to x = 1 with (2.10) leads to 



'1 - 0) min 2 (x) < [ 2 (x)dx 

^'6 [0,1] Jo 



< J J \Vp\ 2 dxdy < C{K + 1). 

[0,1] xT 1 

Let now x* £ [0,1] be any point where 2 (x) attains its minimum on this interval; 
along the particular line x = x* the last inequality yields 



0(x*) < y/C(K+l) (2.12) 

and these oscillations are therefore controlled uniformly in L (C depends only on m 
and Ci). Moreover, x* £ [0, 1] and < p x < C\ control p in average from below and 
from above 

K = Jp(0,y)dy< J p(x*,y)dy<K + c 1 x*<K + c 1 . (2.13) 

T 1 T 1 



For K large enough but fixed (2.12), (2.13) imply < Ki < p(x*,y) < K 2 as desired 



with Ki w K — 0{yK) and K 2 ~ K + (^(vi^) up to constants depending only on 
ci and m. Finally x* £ [0, 1] may depend on L, 5, c\ (and actually does) but stays far 
enough from both boundaries, so that the new translated domain still grows to infinity 
in both directions when L — > +oo. 

□ 
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3 (^-solutions on the infinite cylinder 



From now on we will work in the translated frame Dl =] — L — x*, L — x^xT 1 , where 



x* = x*(L,S) is defined as in proposition 2.3 above. Since the domain depends on L, the 
solution depends on L as well. We emphasize that by writing p = p L (5 > is fixed so we 
may just omit the dependence on 5), and let also set D = K x T 1 to be the infinite cylinder. 

Theorem 3.1. Up to a subsequence we have p L — >■ p in Cf oc (D) when L — > +oo, where 
p G C°° (D) is a classical solution of —mpAp + (c + ct)p x = \ Vp\ 2 . This limit p satisfies 

1. 0<p x < ci 

2. p > 5 

3. p is nontrivial: K\ < p(0,y) < K 2 



where K 1 , K 2 are the pinning constants in proposition 2.3 



Proof. Using interior L q elliptic regularity arguments for fixed q > d = 2 we will obtain W 3,q 
estimates on p L , and this will allow us to retrieve the strong convergence p L — > p in Cf 



loC 



The most difficult term to estimate is |Vp| . We handle it using a different unknown 



which appears very naturally in the original setting (AD-E), namely 



m 2 m+i I m + 1 \ m 

w := p— = m T m+1 (3.1) 

m + 1 \ m J 

An easy computation shows that this new unknown satisfies on Dl a classical Poisson equa- 
tion 

Aw L = f L , (3.2) 

where the non-homogeneous part 

/ L :=(c + a)(/)-^^ (3.3) 

involves only p L and pi", on which we have local L°° control uniformly in L. Indeed, p L is 
pinned at x = by K\ < p L (0, y) < K 2 and cannot grow too fast in the x direction because 
of < pL < c x . 

If m < 1 the exponent — — 1 in (3.3) is positive and we control f L uniformly in L on any 



compact set. However, if m > 1, this exponent is negative and we need to bound pi away 
from zero uniformly in L. For 5 > fixed this is easy since we constructed p L > p~ > 5 > 0, 
but this will be a problem later when taking the limit 5 — > (see next section). 

As a consequence, for any fixed q > d = 2, f L is in L q on any bounded subset Q C D 
and we control 

\\f L \\Li(n) < C 

uniformly in L (C may of course depend on Q, q and 5). Since w L is defined as a positive 
power of p L and p L is locally controlled in the L°° norm uniformly in L the same holds for 



w L , 



\w L \\ L i(u) < C. 
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Let Q =] -a.ajxT 1 C D and K = Q; let also Q 2 =} -2a,2a[xT 1 and Q 3 =} -Sa^afxT 1 
so that Q CC ^2 CC Q$. By interior L q elliptic regularity for strong solutions (the version 
we use here is |14j . Theorem 9.11 p. 235) there exists a constant C depending only on Q2, ^3 
and q such that 

\\w \\w 2 >9 (n 2 ) < C (\\w \\ L i(n 3 ) + \\f L \\Li(n 3 )) ■ 
As discussed above we control w L and f L , hence 

||w L ||vy 2 .9(n 2 ) < C 
for some C > depending only on Q3, Q2 and q. 



(3.4) 



The next step is using (3.1)-(3.3) to express / only in terms of w 



j C + Ol j m 



Expressing V/ L only in terms of w L , Vw L and D 2 w L (which are controlled by | \w L \ \w2, 



using the lower bound p L > 5 > and uniform control on p L , (3.4) implies that 

||V/ |U«(n 2 ) < C 



for some C depending only on the size a of Q. Differentiating (3.2) implies 

A(d lW L ) = d t f L , 



1,2. 



Repeating the previous L q interior regularity argument on CC fi 2 yields 

\\diW L \\ W 2, q{n) < C (||<9^ L || Ln n 2 ) + \ \dif L \\ Lq{ n 2 )) < C, 



and our previous estimate for V/ together with (3.4) finally yield the higher estimate 



\W \\w 3 'i{Q) 



< C. 



The set K = Q = [—a, a] x T 1 is bounded and the exponent q was chosen larger than the 
dimension d = 2. Thus compactness of the Sobolev embedding W 3,q (Q) CC C 2 (K) implies, 
up to a subsequence, that 



L C2 ( K J 

w — > w 



when L — > +00. By the diagonal extraction of a subsequence we can assume that the limit 

'loc 



w does not depend on the compact K. It means w L — > w in Cf oc on the infinite cylinder D. 



The algebraic relation (3.1) and p > 8 > imply that 



J9 > p. 



It implies that we can take the pointwise limit in the nonlinear equation. The limit p solves 
therefore the same equation —mpAp + (c + a)p x = \Vp\ 2 on the infinite cylinder. 

The remaining estimates are easily obtained by taking the limit in < p% < c\, 5 < p~ < 
p L and in the pinning proposition 2.3 Lastly, p is smooth by classical elliptic regularity. □ 
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Proposition 3.1. We have lim p(x,y) = 5 uniformly in y. 

x—>—oo 



Proof. The previous lower barrier 6 < p L on Dl immediately passes to the limit L — > +oo, 
and 

\/(x,y)eD, p>5. (3.5) 
In order to estimate p from above let us go back to the untranslated frame x G [-L, L] 



and remark that by definition p + does not depend on L, see (2.4). An easy computation 



shows that p + (—L) — > 5 when L — > +oo. The subsolution p~ actually depends on L through 
bound? 
that p~ 



boundary condition, see (2.6), but using the monotonicity d x p > is is quite easy to prove 
-L) 



8 when L — » +oo. The left boundary condition consequently reads 



p L (-L,y) = A 



p+(-L)+p-(-L) 



~ 5. 

L— s>+oo 



lim p L (—L,y) is not clear because the convergence 

L— s>+oo 



However, the limit lim p(x, y) 

p L — > p is only local on compact sets (and also because we translated from one frame to 
another) . 

In order to circumvent this technical difficulty we move back to the translated frame and 
build for x e] — L — x*, OfxT 1 a family of planar supersolutions p e (x) independent of L such 
that p £ {— oo) = 5 + e. The construction is the following: fixe > and define p £ (x) as the 
unique solution of Cauchy problem 



p £ (x) : 



—muu" + Cqu' 
u(0) 
u(—oo) 



{u'f 
5 + e 



(3.6) 



where K 2 is the constant in proposition 2.3 such that p (0, y) < K 2 . As already computed 



the setting C = cq < c + a in (3.6) implies that p £ is a supersolution. By monotonicity both 
p L and p £ satisfy condition (2.3), and for L large and S, e small it is easy to check that p <p £ 
on the boundaries x = —L — x*, 0: Theorem 2A_ guarantees that 

V(x,y) E]- L- x*,0[xT\ p L <p £ . 

For S, e fixed, p £ is independent of L: taking the limit L — > +oo yields 

W(x,y) e]-oo,0[xT\ p(x,y) <p £ (x). 



(3.7) 



Taking now the limit e — > in (3.6), it is easy to show that p £ (x) — > p(x) uniformly on 
] — oo, 0], where p is the solution of the same Cauchy Problem as p £ - except for p(— oo) = 5 
instead of p £ {— oo) = 5 + e - and satisfies lim p(x) = 5. Combining the limit e — > in (3.7) 



with the lower barrier (3.5) we finally obtain 



5 < p(x, y) < p(x) 



as desired. 



□ 
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Remark 3. The proof above actually implies a stronger statement than lim p(x, y) = 5, 

x—>—oo 

namely 5 < p < p for x — > —oo: just working on the ODE —mpp" + Cqp' = (p') 2 it is 
straightforward to obtain the exponential decay \p — 5\ = O \p<^ x /' mS \ The exponential rate 
co/m5 degenerates when 5 — > + , which is consistent with the fact that a free boundary 
appears in this limit (see next section). 



As stated in Theorem 3.1 the limit p is non-decreasing in the x direction (as a limit of 
increasing functions p L ). We establish below the strict monotonicity. 

Proposition 3.2. p x > on the infinite cylinder. 



Proof. The argument is very similar to the proof of Proposition 2 J 2_ differentiating the 
equation for p with respect to x yields a linear uniformly elliptic equation satisfied by q = 
Px > 0. The classical Minimum Principle implies that either q > everywhere either q = 0, 
and latter would contradict p(—oo,y) = 5 < K 2 < p(0,y). □ 



4 Limit (5 — > and the free boundary 

In the previous section we constructed for any small 5 > a nontrivial solution p = lim p L 

L->+oo 

of —mpAp + (c + a)p x = \Vp\ 2 on the infinite cylinder D = R x T 1 , satisfying the uniform 
ellipticity condition p > 5 > 0. Let us now write p = p s in order to stress the dependence on 
5. The next step is now to take the limit 5 — > + (S is an elliptic regularization parameter), 
yielding the desired viscosity solution. 

Viscosity solutions are defined as follows: for 5 > let Eg C C°(D) be the set of positive 
smooth solutions p satisfying 



1. lim p(x,y) = 5 uniformly in y. 

x—^—oo 

2. p(x,y) ~ cx uniformly in y at positive infinity. 



According to our Comparison Theorem 2.1 such solutions satisfy p > 5, and let us recall 



that we refer to those as 5-solutions or ^-approximations. We define viscosity solutions as 

Definition 1. A function p G C°(D) is a viscosity solution if there exists a sequence (p S ) s>0 £ 
Ex such that lim p s = p in C?„JD). 

Let us comment on this definition, which is not the standard definition of viscosity 
solutions for second order equations: we will see below that this choice ensures uniqueness 
(see last section), but also regularity. Adapting the previous L q interior regularity argument, 
any viscosity solution p will turn to be C°° on its positive set D + = {p > 0} (see proof of 



Theorem 4.1 below) which is not clear with the usual definition (in addition to being a difficult 
question, see e.g. [7]). This definition can be seen as a particular case of the evanescent 
viscosity method, where we do not regularize the problem by modifying the equation itself 
as it is usually done. The regularization is actually performed through boundary conditions 
p s (—oo,y) = 5 > and p s (x,y) ~ cx, thus ensuring uniform ellipticity p 5 > 5 > 0. The 

+oo 

delicate point is of course the loss of ellipticity when 5 — > + . 
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Our definition [T] includes boundary conditions: any proper setting should consider the 
definition of the notion of solutions independently of boundary conditions. For the sake of 
simplicity we keep this definition, but let us point out that it can be relaxed into a proper 
definition. 

Anticipating that p = lim p s will have a free boundary, we cannot hope convergence to 

hold in any C k topology (k > 1) because of a potential gradient jump across the interface. 
In order to apply Arzela-Ascoli Theorem we need bounds for p s , Vp 5 uniformly in S. At this 
stage we have pinned < Kx < p s (0, y) < K 2 , and < p x < c\ holds on the infinite cylinder: 
we therefore control p s and p s x uniformly on any compact set, but we still have no control at 
all on p s y : 

Proposition 4.1. For any a > there exists C a > such that, for any small 5 > 0, 

x<a^ \p 5 Jx,y)\ < C a . 



Proof. We will first establish this estimate for p L on finite domains [-L — x*,a] x T 1 by 
controlling q = p^ at the boundaries and estimating the value of any potential interior 
extremal point. Taking the limit L — >■ +oo will then yield the desired estimate for p s = 
lim p L . 

L— >+oo 

• Fix a > 0: the uniform pinning K\ < p L (0,y) < K 2 and monotonicity < p x < C\ 
allow us to control p L uniformly in 6, L from above and away from zero on any small 
compact set K — [a — e, a + e] x T 1 . Applying the previous L q interior elliptic regularity 
for w = ^iP* 1 ™ 1 on the slightly larger set ^2 =]o — 2e,a + 2e:[xT 1 DD := K we 
obtain 

llu^llw^n) <C (\\w L \\ Lq{a2) + \\f L \\ Lq{ n 2 )) <C a \\p L \\c HK ) < C a 

for some constant C a depending only on Q, Q 2 and q > 2 fixed, hence on a. It is 
here important that p L is bounded away from zero uniformly in 5 on Q 2 , see proof of 
Theorem 3^ for details. In particular 

\pi(a,y)\<C a (4.1) 

and the monotonicity estimate < p x < C\ combined with the pinning yield 

x < a => <p L {x,y) <C a . (4.2) 



Differentiating (1.4) with respect to y we see that q L '■— Py satisfies the linear elliptic 
equation 

- mp L Aq L + [(c + a)q x - 2Vp L ■ Vq L ] - (mAp L ) q L = -a y p L x . (4.3) 

Let fl a =} — L — x*, afxT 1 : on the left x = — L — x* we had a flat boundary condition 
p L (—x * —L,y) = est so that q L (—L — x*,y) = 0, and on the right boundary x = a 
(4.1) holds. We therefore control \q L \ = Ip^l < C a on the boundaries. 
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In order to control p^ inside fl a we remark that any interior maximum point satisfies 
q > (unless by periodicity p^ = 0, which is impossible if the flow a(y) is nontrivial), 
and of course Aq L < 0, Vg L = 0. At such a maximum point (4.3) immediately yields 

- (mAp L ) q L < -a y p L x - 

using —mp L Ap L = |Vp L | 2 — (c + a)p x as well as the monotonicity estimate 

i 2 

(q L ) 3 - c\q L < \Vp L \ -(c + a)p x q L 

< - (mp L Ap L ) q L < -a y p L p x < C a . 

Since at a maximum point q L > this controls any potential maximum interior point 
max q L (x,y) < C a uniformly in L,5. A similar computation controls q L at any 

potential negative minimum point min q L (x,y) > —C a , and combining with the 

previous boundary estimates yields 

(x,y)e[-L-x*,a]xT 1 => \pfa, y)\ < C a . (4.4) 



Theorem 3.1 ensures that the convergence p L — > p s holds in Cf oc (D): taking the limit L 
+oo in (4.4) finally yields the desired estimate for p 5 . □ 

We can now state the main convergence result when 5 — > + : 

Theorem 4.1. Up to a subsequence we have p 5 — > p in C® oc (D) when 5 — » + , where p > 
is continuous and nontrivial 7^ D + := {p > 0}. Further: 

1. p is Lipschitz on any subdomain ] — 00, a] x T 1 (the Lipschitz constant may depend on 
a). 

2. p solves —mpAp + (c + a)p x = |Vp| 2 in the viscosity sense on the infinite cylinder; 
p\d+ ^ C°°{D + ) is moreover a classical solution on D + . 

3. < p x < ci on D + . 

4- p has a free boundary T := dD + 7^ and there exists an upper semi-continuous function 
I{y) such that p(x, y) > <^ x > I(y) . 

5. If L(yo) '■= liminfJ(y) < I(yo), then at y = y Q the free boundary is a vertical segment 
y^yo 

r n {y = y } = [1(2/0), I(yo)} x {y = y }. 

Proof. The pinning K\ < p s (0,y) < K2 and monotonicity < p s x < c\ control p s on any 
fixed compact set K = J— a , a] x T 1 uniformly in 5. On this compact set p y is moreover 
bounded by proposition 44^ Arzela-Ascoli Theorem guarantees that p 5 — > p uniformly on 
K (up to extraction). Once again by diagonal extraction we can assume that the limit does 
not depend on the compact K, which means local uniform convergence 

5 c ?oA D) 
p — )• p. 

p is nonnegative as a limit of positive functions, and non trivial since for example we had 
pinned < K 1 < p 6 (0, y). 
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1. Proposition 4.1 and monotonicity < p s x < c\ yield C a -Lipschitz estimates on ] — 
oo, a] xT 1 uniformly in 5 for p s : this passes to the C z ° oc limit S — > + , and p is therefore 
Lipschitz on any half cylinder ] — oo, a] x T 1 . 

2. p s G C 2 (D) was a classical solution of —mpAp + (c + = |Vp| 2 on the infinite 
cylinder: according to our definition [l] we need to check that p s grows as cx when 
s — > +oo. This is true but the proof is a long and technical computation: we will prove 
instead that the limit p itself grows linearly (see section [5]). The proof of the linear 
growth is however exactly the same for p s and p, and we therefore admit here that the 
5-solutions grow linearly: p is therefore a viscosity solution in the sense of Definition 

m 

Remark 4. Regardless of this linear growth issue, the limit p is a viscosity solution in 
the classical sense as a consequence of usual stability theorems (see e.g. IT3f §6). This 
is just the classical construction of evanescent viscosity solutions since we had uniform 
ellipticity p 5 > 5 > 0. 

In order to prove the convergence p s p above we could not apply the same local L q 



interior elliptic regularity argument as in the proof of Theorem 3.1 mainly because 

i-^ ninmi ^'i L /ii vi / ii" / 1 ( /-if -i~l"i/-i n/irf(if nrn iii'Vifiiiinit"! ± 

III 



we needed to bound p away from zero (cf. the negative p exponents 1 for the 



non-homogeneous term in (3.2)). This is of course impossible on the whole cylinder 
uniformly in S because the equation degenerates when 5 — > + (this is indeed consistent 
with p = to the left of the free boundary, as claimed in our statement). 

This strategy is however still efficient on the positive set D + = {p > 0}: indeed for 
any fixed compact subset K C D + we know a priori that the limit p is positive, and 
therefore so is p s uniformly in 5 — > + . This allows us to bound p s away from zero 
uniformly in 5 on any compact set K C D + as 

P 5 \ K >C K >0, 

where Ck depends only on K. The interior L q regularity argument in the proof of 



Theorem |3.1| then applies to the letter, and 

p s p in Cf oc (D- 



The limit p G C 2 (D + ) is moreover a classical solution on D + , and smooth by standard 

elliptic regularity. The previous conv 
C°°(D + ) is a classical solution on D~ 



C° (D) 

elliptic regularity. The previous convergence p 5 lo -^ p finally implies that p\o+ = p G 



Convergence p s — > p is strong enough on D + to pass to the limit in < p 5 x < ci, so that 
< p x < c\ on D + . The strict monotonicity is obtained just as for the 5-solutions: 
differentiating the equation for p with respect to x yields and elliptic equation Lq = 
satisfied by q = p x > on D + (where p > is smooth). Applying the Minimum 
Principle shows that either q > 0, either q = 0. Item [4] below will show that p actually 
vanishes far enough to the left: the pinning < K\ < p{0,y) then implies that p has 
to increase at least somewhere in D + , therefore excluding the case q = 0. 



17 



4. In order to show the existence of the free boundary T = d{p > 0} ^ we build 
new suitable planar sub and supersolutions p s ~ (x) , p 5 ' + (x) for p s as follows: defining 
pS,-^s,+ k e un jq Ue planar solutions of the following Cauchy problems 



-muu" + C\u' = (u') ( —muu" + cqu' = {u') 

p 5 '~(x) : { u(-oo) = f , p s ' + (x) : < w(-oo) = 25 

u(0) = K x { w(0) = K 2 

and we have of course $ (p s '~) < $ (p 5 ) = < $ (p 5 ^). Let us moreover recall from 
proposition 3.1 that lim p s (x,y) = 5 uniformly in y, so that p s ~ < p 5 < p s,+ when 



x — > — oo. On the right boundary we set p 6 ' (0) = K\ < p s (0,y) < K 2 = p 5 ' + (0): 
applying Theorem 2.1 on ] — oo, 0] x T 1 yields 



x < => p s '~(x) < p s (x,y) < p s > + (x) (4.5) 



(note that p% ,Px + ,p s x > so that condition 2.3 does hold). When 5 — > + one can 
prove that 



p 5 ' (x) ->■ p (x) := [Ki + Cix] + p s,+ (x) -)• p + (x) := [7T 2 + c x] 



uniformly on IR~, where [.] + denotes the positive part. Taking the limit 5 — > in (4.5) 
yields 

x < =>- p~ (x) < p(x, y) < p + (x) . 

In particular 

x < Xo := =>- p(x, y) < p + (x) = 0, 

c 

K i . . — / \ 

x > xi := =>- p(x, y) > p (x) > 0, 

Cl 

and p has a non-trivial interface of finite width T := <9{p > 0} C {xo < x < xi} as 

pictured in Figure [TJ 

For any y G T 1 the quantity 

7(y) := inf (x G R, p(x, y) > 0) (4.6) 

is well defined because p is nondecreasing in x, and by monotonicity there holds 
p(x,y) > x > /(y). This function /(.) is upper semi-continuous, since its hy- 
pograph 

{(x,y), x < I(y)j = |(x,y), < o| = |(x,y), p(x, y) = o| = D \ D + 

is a closed set (p is continuous). 

5. Assume that y G T 1 is such that L{Uo) '■= lim inf < I(yo)] we prove by double inclusion 

y^yo 

that, if T = d{p > 0}, then V f] {y — y } = [I{y Q ), I(y )} x {y = y }. We write for 
simplicity r := T fl {y = y Q }, and let us point out that by definition p(x, y ) = holds 
for x < /(i/o) • 
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Figure 1: Existence and width of the free boundary. 



r C [I(y ),I(y )] x {y = y } | If x > I(y ) we have that (x Q , y ) G D + , there- 
fore (20,2/0) 4- r = D+/D+ and thus T c] - 00, 1(y )) x {y = y }. If x < l(y ), 
assume that there exists a sequence (x n ,y n ) — > (xo,2/o) such that p(x n ,y n ) G -D + : 
by definition of J(.) we have that p(x n ,y n ) > =>■ x n > I(y n ), and as a conse- 
quence x > liminfJ(y) = /(j/o)- This is impossible since we assumed x < L(yo), 

hence T C [I(y Q ), I(y )} x {y = y }. 



r D [I(y ),I(y )] x {y = y } | Choose any point (x ,y ) G [1(2/0)^(2/0)] x 
{?/ = 2/o} : since p(x ,2/o) = and T = D + /D + we only need to build a sequence 
(x ni y n ) — > (x ,2/o) such that (x n ,y n ) G Let j/ n — >■ y be any sequence such 
that /(2/n) L(Vo)- If ^0 — L(yo) define x n := J(2/ n ) + 1/n: we have that x n > 
J(z/n) =>• p(x n ,y n ) > hence (x n ,y n ) G £>+, and clearly (x n ,y n ) -> (L(y ), 2/o)- If 
now a;o > K2/0), define x n := £o : for n large enough we have again x n > I(y n ) 
hence (x n ,y n ) G D+, and (x n ,y n ) -> (x ,2/o)- Therefore [I(y ), I(y )) x {y = 
2/o} C r . 

□ 



We prove now the Lipschitz regularity stated in Proposition 1.1 



Proof. Under the non- degeneracy hypothesis Px\d+ > a > we prove that the graph of I(y) 
can be obtained as the uniform limit of the £-levelset of p when e — > + , and that these 
levelsets are Lipschitz uniformly in e. 

Let us recall that p\ D + G C°°(D + ): the strict x monotonicity and the Implicit Functions 
Theorem show that, for any e > 0, the £-levelset of p can be globally parametrized as a 
smooth hypersurface 

p(x,y)=e <^> x = I £ (y), 
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Py 

where I £ G C^fT 1 ). Moreover ^ = , and the non-degeneracy hypothesis combined with 

V Px 



proposition |4.1| guarantee that 

dI < <c 



dy 

for some constant C independent of e. By Arzela-Ascoli Theorem we can assume, up to 
extraction, that J e (.) converges to some J(.) uniformly on T 1 . This limit is of course Lipschitz, 



and we show below that J(y) = I(y), where I is defined as in Theorem 4.1 (p(x,y) > <^ 
x > I{y)). 

By continuity we have that p(I e (y),y) = e =>- p(J(y),y) = 0. Chooss xq < J(y) and e 
small enough: integrating p x > a > from x = I e (y) < Xo to x = xq leads to p(xo,y) > 
e + a(xo — I e {y))- Taking the limit e — > + yields p(xo, y) > a{xo — J(y)) > 0, and therefore 



p(J(y),y) = 
x > J(y)=> p(x ,y) > 



| J(y) = inf(x, p(x, y) > = I(y) 



by definition (4.6) of /. Thus I = J is Lipschitz, and by continuity T = {x = I(y)}. □ 



Proposition 4.2. The corresponding temperature variable v = € C(D) solves the 

original equation A(t> m+1 ) = (c + a)v x in the weak sense on the infinite cylinder D: for any 
test function \l/ e T>(D) with compact support K C D we have that 

v m+1 /\m<lx<\y + J J (c + a)v^ x dxdy = 0. 

K K 

Proof. We denote by v L and v 6 the temperature variable corresponding to our two successive 
approximations p L and p s . Let $ G D be any such test function with compact support 
K C D: let us recall that the finite cylinder grows in both direction, and consequently 
K C Dl for L large enough. p L > was a smooth solution of —mpAp + (c + a)p x = \Vp\ 2 
so that v L was a smooth solution of A(v m+1 ) — (c + a)v x = 0, and therefore 

(v L ) m+l A^dxdy + ^ (c + o> L ^daxh/ = 0. 

K K 

When L — > +oo the Cf oc (D) convergence p L — >■ implies the C° oc (D) convergence f L — >■ f 5 , 
hence 

y ^ (u 5 ) m+1 A^dxdy + J J (v s )(c + a)w^ x dxd?/ = 0. 

K K 

Using the C z ° oc (D) convergence p 5 — > p the integrals above finally pass to the limit (5—7-0. □ 
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5 Behavior at infinity 



We prove in this section that the behavior at infinity is not perturbed by the shear flow, 
compared to the classical PME traveling wave p(x,y) = c[x — xq] + . As mentioned above 
the results of this section are established directly for the final viscosity solution p = limp 5 , 
but easily extend to the (^-solutions. 

Theorem 5.1. p(x, y) is planar and x-linear at infinity, with slope exactly equal to the speed: 

p x (x, y) ~ c p y (x, y) ->■ 0, p(x, y) ~ cx 
uniformly in y when x — > +oo. 

We start by showing that p(x, y) grows at least and at most linearly for two different 
slopes; using a Lipschitz scaling under which the equation is invariant, we will deduce that 
p is exactly linear and that its slope is given by its speed c > 0. This will be done by 
proving that in the limit of an infinite zoom-out (x, y) — » (X, Y) the scaled solution P(X, Y) 
converges to a weak solution the usual PME (a = 0) which has a flat free boundary X = 
and is in-between two hyperplanes. By uniqueness for such weak solutions of the PME our 
solution will agree with the classical planar traveling wave P(X, Y) = [cX] + , hence the slope 
for p(x, y) at infinity. 

5.1 Minimal growth 

Since p x < c\ we have an upper bound at infinity p < Cix; we show in this section that we 
also have a similar lower bound: 

Theorem 5.2. There exists C>0 such that 

x > =^ p(x, y) > Cx. 



Let us recall that we have pinned 

K x < p(0, y) < K 2 , K < I p(0,y)dy<K + C 



T 1 



where K\ > K — CyK and K<i < K + C\/K. The constants C above depend only on m > 
and the upper bound for the flow c\ > c + a(y), and K > can be chosen as large as required 



(see proof of proposition 2.3 for details). 
We will denote by 

0(x) = maxp(x,y) — minp(x,y) 
the oscillations in the y direction, which is a relevant quantity that we will need to control. 
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Lemma 5.1. There exists a constant C > and a sequence {x n ) n>0 G [n,n + 1] such that 



0(x n ) < C 



i 



p(n + l,y)dy 



T 1 



Proof. Integrating by parts —mpAp + (c + a)p x = \Vp\ 2 over K n — [n, n + 1] x T 1 we obtain 
(m — 1) J J | Vp\ 2 dxdy + m J pp x (n, y)dy — m J pp x {n + 1, y)dy + J J (c + a)p x dxdy = 0. 

(5.1) 

We distinguish again m < 1 and m > 1: 



1. If m < 1 we use pp x (n + 1, y) > 0, < c + a < ci and < < c\ in (5.1) to obtain 



(1 — m) J J | Vp\ 2 dxdy < mj pp x (n, y)dy + J J (c + a)p x da;d?/ < mci / p(ra, y)d?/ + c\. 

Choosing if large enough we can assume by monotonicity that c\ < mc\ f p(n,y)dy, 
and therefore 

\Vp\ 2 dxdy < / p(n,y)dy. 

1 — m J 

If x n G [n,n + 1] is any point where 0(x) attains its minimum on this interval, then 

n+l 

0(x n ) < J 0(x)dx 

n 

n+l 



< I [f\Py(x,y)\dy\dx < C JJ\Vp\ 2 dy. 

n XT 1 / V K n 

Using our previous estimate and monotonicity we finally obtain 



0{x n ) < C 



i 



p(n,y)dy < C 



T 1 



p(n + l,y)dy, 



T 1 



where C depends only on m and c\. 



2. If m > 1 we use pp x (n, y) > 0, (c + a)p x > and p x < c\ in (5.1 ), yielding 



(m — 1) J J \Vp\ 2 dxdy <mj pp x {n + l,y)dy < mci ^ p(n + l,y)dy. 

ii"n T 1 T 1 

The rest of the computation is similar to the case m < 1. 



□ 
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Corollary 5.1. There exists C > such that 



x > 



0(x) < C 



\ T 1 



p(x,y)dy. 



Proof. Since < p x < C\ the function 0(x) is clearly 2ci-Lipschitz, and by Lemma 5.1 
ensures that 



0(x) < 0(x n ) + 2 Cl <C 



\ 



p(n + l,y)dy + 2c x 



for any x G [n, n+1]. By monotonicity we can moreover assume that 2ci < Cy f p(n + l,y)dy 
if K is chosen large enough, and therefore 



0(x) < C 



p(n + 



For the same reason we can also assume that 



p{n + 1, y)dy < / p(x, y)dy + c 1 < C / p(x, y)dy 



and combining with the previous inequality yields the desired result. 
Proposition 5.1. For any x > we have that 



d I f m+1. \ 771 + 1 

dx [ Pm ^ V)dy , ,„ 



[c + a{y))p m {x,y)dy. 



T 1 



□ 



Proof. We establish this equality for the uniformly elliptic solution p s > 5 up to a constant 
Cg, with Cg — > when 5 — > 0. The equation for p s can be written in the divergence form 

V- (V)™ V/) = ((c + a) (p 5 )™) , 

and integrating by parts over = [xi,X2] x T 1 yields 

(p 5 ) ™ p*(a? 2 , y)rfy-y (/) ™ pi(xuy)dy = J (c+a) (p s ) ™ (x 2 , y)dy-y (c+a) (/) ™ (x x , 2/)tfy 

T 1 T 1 T 1 T 1 

for any xi < x 2 . As a consequence, the quantity 



F(x) := j (p S ) ™ p£(x, - j (c + a) (/) ™ (x, = C 5 



(5.2) 
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is constant. Let us recall from proposition 3.1 that p s (—oo,y) = 5 uniformly in y, and also 
the uniform bounds Co < c + a < c\ and < p s x < C\\ taking the limit x — > — oo in (5.2) 

leads to C s = O 1 8 



Fix any x > 0: the strong C\ oc convergence p s — > p on D + = {p > 0} is strong enough to 



take the limit 8 — > in (5.2), which reads 



P m Px(x, y)dy - j (c + a)pm (x, y)dy = 

T 1 T 1 

Finally, p is is smooth for x > (because p > 0), and the last equality above easily yields 
the desired differential equation. □ 



We can now prove the claimed minimal growth: 

m + l 



Proof, (of theorem 



5.2). Define f(x) := J p ™ ( y x,y)dy: Proposition 



T 1 



f'(x) = m + 1 j (c + a)p™dy 
m J 

T i 



5.1 



reads 



for x > 0. By monotonicity j p(x,y)dy > j p(0,y)dy = K, and by corollary 

T 1 T 1 



5.1 



(5.3) 



we control 



the oscillations 0(x) < Cy f p(x,y)dy. Choosing K large enough the oscillations of p are 
small compared to its average along any line x = est > 0, hence 



m+l 



(c + a)p m dy > c J p m dy > C | / p ™ dy 

T 1 T 1 XT 1 



Cfm+i (x) 



This estimate combined with (5.3) leads to f'(x) > Cf m + l (x), and integration yields 

f™+*(x) > Cx. 
Finally, since we control the oscillations of p, 



p(x,y) > C J Pi.x } y)dy >C\Jp™ (x,y)dy | > Cf™+i(x) > Cx. 



□ 



5.2 Proof of Theorem 5.1 



We start by estimating how fast p becomes planar at infinity: 
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Proposition 5.2. Let as before 0(x) := max p(x,y) — min p(x,y); there exists C > such 
that when x — > +00 

0(x) < -. 

x 

Proof. For x large enough we know that p(x,y) > is smooth; w := -^ip"^ is therefore 
smooth, and satisfies as before 

= /, / = (c + a)p™ -1 ps- 

We will first show that the y oscillations of w cannot blow too fast when x — > +00, and then 
deduce the desired planar behavior for p. 
The Fourier series 

w(x,y) = J2Mx)e 2lnny 
is at least pointwise convergent, and for n^Owe have that 

-<(*) +i7r 2 n 2 w n (x) = f n (x), f n (x) :=-J f{x,y)e- 2 ^dy. (5.4) 

T 1 

The oscillations of w in the y direction are completely described by its Fourier coefficients 



w n (x) for n 7^ 0, in which case (5.4) is strongly coercive. This coercivity will allow us to 
control how fast w n (x) may grow when x — » +00, and therefore how much w can oscillate. 
Since p is at least and at most linear and p x , c + a are bounded we control 

\f n \(x) ^Cx^ 1 (5.5) 



uniformly in n. Moreover, taking real and imaginary parts of (5.4), we may assume that 
w n (x), f n (x) are real and that n = \n\> 0. 

• We claim that there exists C > such that, for any n^O and x — > +00, there holds 

C 1 

\w n (x)\ < —x™~ x . (5.6) 



7?. 2 



Indeed, since < w = JB ^rv'"'^ < Cx"m we have that 

|2/„\ / IL../„ MI2 / r<V- 



\w n \ (x) < \\w(x, OIIz^t 1 ) <Cx z 



As a consequence w n cannot have a component on the homogeneous solution e +27rnx of 



(5.4) for and it is then easy to see that it is explicitly given by 



w n (x) = e- 2 ™ ix - xo) w n (x ) + e~ 2nnx j e^ nz | / e~ 27rnt f n (t)dt | dz. (5.7) 



Our claim (5.6) is then easily obtained manipulating this explicit formula, the com- 
putations involving several integrations by parts and the fact that w n (xo) is rapidly 
decreasing in n (since w(xo, .) G C°°(T 1 )). 
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As a consequence of (5.6), the series 

J w(x,y)dy = Y,Mx)e 2i7vny 



w (x,y) := w(x,y) - 



ji 



is uniformly convergent and |w _L (a;,?/)| < Cx™ 1 . This clearly bounds the oscillations 
of w when x — > +00 by 

max w(x,y) — min w(x,y) < 2\\w ± (x, OIU^fT 1 ) < Cx™~ x . (5.8) 
Translating the oscillations of w in terms of those of p = Ciu™^ 1 leads to 



0(x) = C x ( max u>™+i (x, y) — min tym+i (x, y) 

j/ST 1 ysT 1 



< C ( min w(x, y) 



m 1 

m + 1 



max w(x,y) — min w(x,y) 



Since w = ^jp'"" 1 > Cx™™ 1 and — 1 = — estimate (5.8) finally implies that 

1 



0(x)<c{x m ^y m+1 xCx^ 



x 

□ 



For any e > let us introduce the Lipschitz scaling 



P £ (X,Y)=ep(x,y), (x, y) = F); 

when £ — ?• + this corresponds to zooming out on the whole picture. Uppercase letters will 
denote below the "fast" variables and functions, whereas lowercase will denote the "slow" 
ones. Since we want to zoom out it will be more convenient to consider below the cylinder 
D = M x T 1 as a plane M 2 with a 1-periodicity condition for p in the y direction, correspond- 
ing to a plane with e-periodicity in Y for P £ . 



The proof of Theorem 5^ relies on three key points: the first one is that the equation is 
invariant under this scaling. The second one is that, since the shear flow a(y) is 1-periodic 
with mean zero, the corresponding flow A e (Y) = a(Y/e) is e-periodic with mean zero in 
Y: Riemann-Lebesgue Theorem guarantees that A £ — ^ in a weak sense when e — > 0, so 
that any limiting profile P = limP e will not "see the flow" and thus satisfy the usual PME 
—mPAP + (c + 0)Px = I VP| 2 . Finally, proposition 5.2 guarantees that the oscillations of 
p in the y direction decrease at infinity: zooming out, the limit P will therefore be planar, 
P Y = 0. 

In the limit of this infinite zoom-out the scaled profile indeed converges: 

Proposition 5.3. Up to a subsequence we have P e (X,Y) — > P(X,Y) when e — > + . The 
convergence is uniform on IT x E and C} oc on IR + * x M.. Further: 
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1. P is continuous on the whole plane and P = for X < 



2. < CX < P(X, Y") < ciX for X > 0, where C > is the constant in Theorem\5.2 
and C\ > c + a(y) is the upper bound for the flow. 

Proof. We pinned the original solution p such that < p(x,y) < K 2 for x < 0, and this 
immediately implies that P £ = ep < eK 2 — > uniformly on the closed left half-plane X < 0. 
On the right half-plane < P|-(X, Y) = p x (x, y) < Ci bounds P £ from above as 

P £ (X, Y) < P £ (0, Y) + Cl X < K 2 e + c x X, (5.9) 



and Theorem 5^ bounds P £ away from zero 

P £ (X, Y) = ep(X/e, Y/e) > CX. (5.10) 

Let us recall that p is a smooth classical solution on D + = {p > 0} D 1R + xT 1 : for £ > 
the rescaled profile P £ is therefore a smooth classical solution of the rescaled equation 

-mP £ A xx P £ + [c + A £ (Y)] P £ x = \V x ,yP £ \ 2 , A £ (Y) = a{Y/e), 

at least for X > 0. 

Using our previous interior elliptic L q regularity argument for 

2 

W e := jn F £ : = ( c + v4 £ )(P £ )-" 1 Pf : , A1U £ = P £ , 

771+1 

we obtain as before an estimate 

I 1^11^2.9 (Boo) < C 

uniformly in e on any ball B\ C M + * x R of radius 1 and for q > d = 2 (see proof of Theorem 



3.1 for details). It is here important that P £ is bounded away from zero uniformly in e for 



X > 0, see again proof of Theorem 3.1 (in particular the case m < 1). By compactness 
W 2,q CC C 1 on bounded balls (q > d = 2) and moving the center of the ball B\ we may 
assume, up to extraction of a subsequence, that 

W £ — > W in C} oc (R + * x R). 

Since we took care to step out of the zero set uniformly in e, this convergence easily translates 
into 

Cj (R+*xR) 



and P is continuous on IR + * x M as a locally uniform limit of continuous functions. Taking 
the limit e ->■ for X > in CX < P £ (X, F) < K 2 e + c x X we obtain 

X > CX < P(X, Y) < ciX 

as claimed, which gives as a by product the continuity along X = (let us recall that P = 
on the left half-plane). □ 
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Remark 5. No higher regularity can be obtained with this interior elliptic regularity argu- 
ment: C 2 convergence would require for example W 3,q estimates, involving V ' ( X y)F £ which 
contains the singular derivative dyA £ = \d y a. 

As usual we need to determine the limiting equation satisfied by the limiting profile in 
some sense: 

Proposition 5.4. The limiting junction P solves the PME 

-mPA {XtY) P + cP x = |V (x ,y)P| 2 
in the weak sense on the whole plane. 

Proof. By definition of weak solutions we want to prove that, for any test function $(X, Y) 

with compact support K C M. 2 , the corresponding temperature V(X, Y) := (^j-j-P(X, Y)j m 
satisfies 

I := J J V m+1 A$dXdY + J J cV<$> x dXdY = 

K K 

(note that the shear flow A £ (Y) <-> a(y) disappeared in the advection term). Let us recall 



from Proposition 4.2 that p was a weak solution on the whole plane, and that the equation is 



invariant under Lipschitz scaling: for any e > the scaled temperature V e therefore satisfies 
1(e) := J J (V £ ) m+1 A$dAW + Jj(c + A £ )V £ <S> x dXdY = 0; (5.11) 

K K 

the problem is as usual to take the limit in this formulation. 

• If K C Mr* x K this limit is straightforw ard: (c + A 6 ) is uniformly bounded (cq < 
c + A £ < Ci), and according to Proposition 5.3 V s = (^^P 6 ) m — > uniformly on K. 



If K C 1R + * x M the limit V is positive so there is no such trivial convergence; it is 



convenient to split (5.11 ) in three parts I — I\ + I\ + 13 — 0, with 

I x (e) := [[ (V £ ) m+1 A$dXdF , 



h(e) := c V e $ x dXdY, 



K 

:= cj ' J V Sfx^ 
K 

J 3 (e) : = J J A £ V £ $ x dXdY. 

K 

'0 



The C[ oc convergence P £ — > P shows that I\ and 1 2 immediately pass to the limit. To 
deal with ^3 we compute with Fubini Theorem 

I 3 (e) =\l A £ V £ $ x dXdY = [ A £ (Y) I / V £ {X, Y)<S> X {X, Y)dX I dY 



K 
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since $ has compact support and V s — > V uniformly on K we deduce that ^> £ {Y) — > 
^(Y) uniformly on M. and \l/ have both compact support: the convergence \l/ e — > ^ 
therefore also holds in L X (]R), and by Riemann-Lebesgue Theorem A 6 — ^ weakly in 
L 1 (IR) (let us recall that A £ (Y) is e periodic with mean zero), /3(e) is therefore a 
dual pairing 13(e) = (A £ ,^ £ )( L t iiL1 ) of a weakly converging sequence with a strongly 
convergent one: hence the limit /3(e) — > 0. 

If K n {X = 0} 7^ the convergence is more delicate because K crosses the free 
boundary and we do not have uniform convergence V s — > V on K\ however since 
P(0, Y) = both V and V s have to be small on a neighborhood of K n {X = 0}. 
For small r > we prove that there exists Eq > such that for all e < Eq there holds 
\I-I(e)\ < r. 

For 1] > to be chosen later let us define the partition 

k = n {X < u (x n {\X\ < rj}) u (x n {X > +77}) ; 

V V V 

:=A-- :=2C! :=^+ 

X' 7 is a striped ^-neighborhood of X fl {X = 0}. On we already proved that 
Ji,/2,/3 converge: we only have to cope with the contribution from K v , and it is 
clearly enough to prove separately 



(v £ ) m+1 - v m+1 \ .|A$|dXdr < 



3 



A'') 



c \V*-V\.\* x \dXdY < § . (5 . 12) 



K'> 

A £ \.\V £ -V\.\$x\dXdY < 



3 



A' 1 ' 



Let us recall the previous bounds for the pressure variables, derived from the scaling 
and the Lipschitz estimate in the X direction: 



-r]<X<0 

0<X<T] 



0<P £ < K 2 e 
P = 

< P £ (X, Y) < P £ (0, Y) + ciX < K 2 e + ciX 
P < ciX 



Choosing 77 and e small, any positive power of the pressures P £ , P can clearly be made 
as small as required on K v ; this is also true for any positive power of the corresponding 
temperatures V s , V (being themselves positive powers of the pressure), and all the 
terms |A$|, \ A £ \ are bounded uniformly in e: we complete the proof using the 



celebrated triangular inequality in the integrals (5.12). 

□ 



We can now finally prove Theorem 5.1 
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Proof. By Proposition 5.3 and up to extraction we have that P £ — > P uniformly on R~ x R 
and locally in C 1 (IR + * x R); the corresponding temperature V > is a weak solutions of 
the stationary PME — Ax,y (V m+1 ) + cVx = (previous proposition) and has a flat free 
boundary X = separating P = to the left from P > to the right, where it is in-between 



two hyperplanes CX < P(X, Y) < C\X. Moreover, proposition 5.2 shows that the limiting 
profile is planar, dyP = 0. Indeed, for fixed X > and any Yi,Y 2 , we have for e small 
enough 

iP^Xo^) - P%X ,Y 2 )\ = e\p(X /e,Y 1 /e)-p(X /e,Y 2 /e)\ 

e 



< e x C 



X ™ 



taking the limit e ->■ yields \P(X , Y 2 ) - P(X , Y 2 ) \ = lim|P £ (X , Y 2 ) - P £ (X , Y 2 )\ = for 

any Xo > and Yi,Y 2 . It is well known that there exists only one such planar solution of 
the PME, which is the standard planar traveling wave 

P(X,Y) = [cX} + . 

Since the limit is unique the whole sequence actually converges, lim P £ = P: for any 

x £ = \ -> +oo the Cf oc convergence P £ (X, Y) -»■ [cX] + on R + * x R shows that 

1 



max \p (x e , y) — cx £ \ = max 

j/GT 1 " Ye[0,e] 



-P £ (ex e ,Y) - cx £ 



= x £ max |P e (l,y)-P(l,y)| = o(x e ), 

which means precisely y) ~ ca; uniformly in y when x — > +oo. Using now the stronger 
C} oc convergence for X > we finally obtain 

max \p x (x 6 , y) - c| = max |P|(l,y) - P x (l,^)| = o (1) =^ p x ~ c 

yeT 1 Ke[0,e] e^O 

max \p v (x £ ,y) - 01 = max IP5 (1, K) - Py(l, F)| = o (1) ^> p„ ->• 0. 
veT 1 yefo,el ' e^O 



□ 



5.3 Asymptotic expansion at infinity 

We have shown that p(x, y) ~ cx uniformly in y when x — > +oo. In this Section we strengthen 
this estimate and derive and asymptotic expansion 

p(x,y) =cx + q{x,y) 

with W 1)0 ° estimates on q as x — > +oo. 
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For any function f(x, y) periodic in the y direction, we denote the average (the projection 
onto constants in L 2 (T 1 )) by 

(f)( x ) : = / f{x,y)dy- 



T 1 

The orthogonal projection onto functions with mean zero is denoted by 

f ± (x,y):=f(x,y)-(f)(x). 

The x derivative commutes with both these projectors, 4:(f) = (fx) and (f x ) ± = (f^x- 
The ansatz p(x, y) = cx + q(x, y) gives 

(p)(x) = cx+ (q)(x), p L {x,y) = q^i.x.y), 

and q, (q), q x are o(x). The main result of this section is 

Theorem 5.3. When x — > +00, we have that: 

1. For any m^l the correction q{x,y) becomes planar: there exists C > such that 

C 



\q ± \(x,y) + \Vq ± \(x,y) < 



x 



2. Assume in addition that 1 < m ^ N* ; and let N = [m]: there exists a finite sequence 
qi, qN £ K and some q* G M. such that 

q(x, y) = x (q\x~™ + q 2 x~™ + ... + qNX~™^j + q* + o(l). 
The orthogonal projection p- L (x,y) is controlled by the oscillations in the y direction 



\p (x,y)\ < O(x) = max p(x,y) — min p(x,y), and Proposition 



5.2 



therefore implies that 



\q ± \(x,y) = \p ± \(x,y)<- (5.13) 

x 

when x — > +00. 

We prove the first estimate of the Theorem as a separate Proposition. 
Proposition 5.5. There exists C > such that 

C 

\q^{x,y)\ + \Wq ± (x,y)\ < -. 

x 



Let us stress that this statement holds for any m, although we will specifically consider 
m > 1 in the sequel. 
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Proof. By (5.13) we already control \q x \, and it is enough to control its gradient. The 

{c + a)p x \Vp\ 2 



equation for p reads 



Ap 



(5.14) 



mp mp 

and when x — > +00 we know that Vp — > (c, 0) and p ~ cx uniformly in y: as a consequence 



< — . Averaging in y yields < — , and therefore 

|A(^)| = |A(p^)[ = [A^ - <p>"[ < 



C 



X 



Choose now xo large and 1/0 G T 1 , and denote by E\ the ball of radius 1 centered at (xo, yo) 
As discussed above there exists C > such that, if xq is chosen large enough, 



{x,y) e B x 



\Aq x \(x,y) < 



c_ 

xo c 



The constants above depend on the radius of the ball R = 1 but not on its center. Finally, 
the classical elliptic theory for Poisson equation on a ball controls the gradient at the center 
by |Vg _L |(xo, yo) < C (\\q ± \\l x '{b 1 ) + I |^Q'" L | |z°°(Bi)) > with C depending only on the radius of 
the ball. □ 

As a corollary, we have that 

Lemma 5.2. If m > 1, there exists A G R such that 

A 



(cx + (q)y< 



- + o 



x- 



(5.15) 



holds when x — > +00. 



This technical result will later allow us to establish the asymptotic expansion q = x( 
stated in Theorem 15.31 



Proof. Equation(5.14) with p(x,y) = cx + q(x,y) leads to 

|Vg| 2 



mAq 



[a - c)q x 



+ 



ca 



By proposition 5.5 we control \q~ 



cx + q cx + q cx + q 
■ O (1/x), and it is easy to expand 



(5.16) 



1 



1 



1 



cx + q cx + (q) + q 1 - cx + (q) 



1 



cx + {q) 



O 



1 



This expansion allows us to estimate separately the three terms in the right-hand side of 
(5.16), and in particular their average in y. 
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The first one is 



+ 



cx+(q) 



(q)' + 



<?>' 



cz+(<z> (cx+{q)f 

C , , , / 1 



-a — 



for + ^ 



cx + (q) cx + (q) \x 3 



purely orthogonal lower order 



Averaging in y then yields 

(A)(x) = 



' (q)'+ Y 



cx+(q) cx + (q) (cx + (q)) 



x- 



(5.17) 



We expand the second one as 



cx+q I 



cx+(q) 



((9)') 2 +|V9 ± r + 2( 9 >'«)J 



purely orthogonal lower order 



and averaging leads to 



(B)(x) 



The last term is 



C(x, y) := —^—ca 



(aq ± )+ C a +0[- 
v ' cx + (?) \ x a 



(5.1* 



(cz+(g»< 



purely orthogonal lower order 



and finally 



(C)(x) 



ziaq 1 ) + ' 1 



(cx + (?)) 



X" 



(5.19) 



Averaging (5.16) in y reads m(q)"(x) = (A)(x) — (B}(x) + (C)(x): taking advantage 



of (5.17)-(5.18)-(5.19) and rearranging, we obtain 



cx + (q) \cx + (?) / V -r 

Multiplying by the integrating factor (cx + (?)) m yields 



m \cx + (?) 



(5.20) 
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If f(x) 



(cx+{q))m f (gqr x ) N ' 
cx+(q) 



denotes the first term in the right-hand side above, an 



integration by parts combined with (g- 1 ! < C/x =>■ \(aq ± )\ < C/x allows us to show 
that / is integrable at infinity and that 



+00 

J f{z)dz = o[x^- 2 ^ 



This is precisely where we used the technical assumption m > 1: otherwise this term 
may not be integrable at infinity. 



Equation (5.20) can therefore be integrated from x to +00: there is a A 6 1 such that 

+00 

{cx + (q))™ (q)' - A = - J f(z)+o(z&- 3> ) dz = o(x^- 2 y 

x 

and we conclude the proof dividing by (cx + (q)) m ~ Cx™. 

□ 



We finally prove Theorem 5.3 



Proof. The first item is stated in proposition |5.5 Regarding the second item, let us recall 
that q = (q) + q ± and that |g" L | + |Vg _L | < C/x: our statement is actually that the asymptotic 
expansion holds for (q) instead of q, since the transversal part {q^l is negligible when x — > 
+00. 



Let us recall from Lemma 5.2 that (q}(x) satisfies 

A 



(cx + (q)Y 



- + o 



x- 1 



(5.2i; 



for some A £ R. If A = then (q)' is integrable and our statement immediately holds with 
qi = ... = q N = 0. 



If A 7^ (5.21) with cx + (q) ~ cx yields (q)' ~ Ai/x™, which is not integrable if m > 1: 



integrating therefore yields (q) ~ qix 1 ™. Injecting this equivalent into (5.21) leads to 



(?>' 



A 



O 



cx + qix m + 



O I X 1 m 



Expanding the quotient in Taylor series at order two in powers of x m yields now 

1 



(<?)' = m + ^2% m + O ( X 



+ 



and integrating 



(q) = x {^qix ™ + q 2 x ™ \ + o {^x 
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Injecting again into (5.21) yields the next order, and so forth: by induction one shows that 



(q) = x (q x x m + ... + q k _i% k ™\ + o (x 1 V j 

(5.22) 

(q)' = \\x~™ + ... + \ k x~™ + o \x~™} + O (x~ 2 ) . 

k 

As long as k < N = [m] < m the last term \kX~™ in the expansion of (q) above is 
not integrable, and we may continue the induction 



(q)' = X\X m +\kX ™+o(x m j+0(x 2 ) =>- (q) = x (q\X m + ... + q^x m j 



+o\x 1 



• If now k = N + 1 = [m] + 1 > m, the terms X k x m + o (^x m j + O (x 2 ) in (5.22) are 
integrable: integrating one last time we obtain as desired 

(q) = x (qix~™ + ... + q N x~™j + q* + o(l), 

where q* is the constant of integration. 

□ 

Remark 6. Let us stress that the condition m N is purely technical. If m = [m] = N is 
integer we may obtain at some point (q)' = \ix~™ +... + — + .. in the induction above. This 
would yield of course a logarithmic term, which would have to be properly taken into account. 
An asymptotic expansion could be obtained nonetheless, but the resulting computations would 
be long and not very insightful. 



6 Uniqueness 

In this section we prove that the wave profiles of <5-solutions are unique for given 5 > (and 
of course up to x-translations). Since we defined viscosity solutions as limits of 5-solutions 
when 5 + , uniqueness of viscosity solutions should follow. In order to keep this paper 
in a reasonable length we will not take this limit, which would require significant amount 
of technical work: one should indeed build a family of (^-solutions (p 5 ) s>0 such that the 
whole sequence converges to a non-trivial viscosity solution p when 5 — > + (so far we only 
extracted a subsequence). 

Let us point out that all the results in Section [5] are stated for the final viscosity solution 
p = lim p s , but easily extend to the ^-solutions for 5 > 0. Through this whole section we 
fix S and denote by p,pi,P2 any (smooth) 5-solutions in order to keep our notations light. 
For the sake of simplicity we only consider the case 1 < m ^ N, for which the asymptotic 
expansion at infinity p = cx + gix 1- ™ + ... holds (see remark ^ above). 

The main result of this section is 
Theorem 6.1. The 5-solutions are unique up to finite x -translation. 
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Let us start with some technical statements: 
Proposition 6.1. Any d~ -solution has an asymptotic expansion 

p(x, y) = cx + x (q\X~™ + ... + q N x~™^j + q* + o(l) 

uniformly in y when x — > +00, where q±...qN, q* £ K and 1 < N = [in] < m. 

Remark 7. The coefficients qi, q* and the remainder o(l) above may depend of course on 
5, which is fixed here. 

Proof. We may proceed exactly as we did for the final viscosity solution, see section [5] and 
in particular the proof of Theorem |5.3 □ 



The following holds at negative infinity, where we recall that p(—oo,y) = 5 > uniformly 
in y. 

Lemma 6.1. We have that 

|Vp| ->■ 0, \D 2 p\ -> 

uniformly in y when x — > —00. 

Proof. Let w := ^jp 2 ™^ and / := (c + recall that the Poisson equation 

Aw = f 

holds in the whole cylinder. Taking advantage of p(— 00, y) = 5 > (w thus being uniformly 
bounded from above and away from zero) we may safely apply our previous interior elliptic 
regularity argument on Q n :—] — n, —n + lfxT 1 (n G N) to show that 



< C 



for some constant C > and q > d = 2 both independent of n. 
Setting 

OHO^fxT 1 , p n (x,y):=p(x-n,y), 
the previous estimate reads 

||p n ||iy3^(n) < C. 

By compactness W 3 ' q (Q) CC C 2 (Q) we may extract a subsequence p Uk — > p°° in C 2 (Q). 
Since p(— 00, = 5, the limit p°°(a:, y) = p(— 00, y) = est = 5 is unique: standard separation 
arguments show that the whole sequence converges 

p n 5 in C 2 (H). 

which immediately implies our statement. □ 
Proposition 6.2. The coefficients q\...qN in the asymptotic expansion are unique. 
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Proof. Let pi and p 2 be two different (^-solutions, thus satisfying 

Pi = cx + x (qi,ix~™ + ... + qi t N%~™ ) + q* + o(l) 

p 2 = CX + X ^g 2 ,l^~™ + ••• + <?2,iV^™) + <? 2 + 

when x — > +00 for some constants q^k, q* £l, £ = 1, 2, = l...iV and A/" = [m]. 

Assume by contradiction that q^i > q^y. we will first slide p 2 far enough to the right so 
that P2 < Pi on the whole cylinder. Slowly sliding p 2 back to the left we will obtain a contact 
point between p\ and a translate of p 2 , thus contradiction the classical Maximum Principle. 



Lemma 6.1 allows us to pin p\ such that, for x < 0, there holds 

1. 6<p 1 (x,y)<pi{0,y)<26, 

2. |Vpi| and |Api| are small. 

This can be done suitably sliding, since p± — > 5, |Vpi| —> and \ Api\ — > when x — > —00. 
In this proof p\ will be fixed once and for all, and we will only slide p 2 with respect to p\ 
(for the sake of clarity p 2 denotes below any translation). 

• Since we assumed that q\^\ > g 2j i we have 

[Pi - Pa] (+00, j/) = +00 

for any (finite) translation p 2 . Using d x pi > we may therefore slide p 2 far enough to 
the right so that 

x > Pi{x,y) > p 2 (x,y). 

We claim that, applying a suitable comparison principle, we may assume that p\ > p 2 
also holds for x < 0. In order to see this, define z '■ = P\—pi and subtract the equation 
for p 2 from the equation for p\ to obtain 

C[z] := —mp 2 /S.z + (c + a)z x — (Vpi + Vp 2 ) • Vz — (mApi)z = 0. (6.1) 

Testing z(x) := e Xx as a supersolution for some A > 0, an elementary computation 
leads to 

C[z] = e Xx (-mp 2 \ 2 + (c + a)X - {d^ + d x p 2 )\ - mA Pl ) . 

Sliding p 2 far enough to the right we have, for x < 0, that p 2 ~ 5 and that d x p2 is 
negligible. Since we also pinned |Vpi| and |Api| to be small, the main contribution in 
the parenthesis of the right-hand side above comes from the first two terms. Choosing 
A > small enough, it is clearly possible to satisfy 

— mp 2 A 2 + (c + a)X <; — m5X 2 + c A > 

(choose for example A = co/2m5), and therefore 

x < => C[z] > 0. 
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Setting z := wz, the new variable w satisfies this time an elliptic equation 



C[w] = 



where £ is uniformly elliptic, has positive zero-th order coefficient C[z] > 0, and 
therefore satisfies the Minimum Principle. 

One the right boundary x = we may assume that Pi(0, y) > P2(0,y) (once again 
sliding p2 far enough to the right), and therefore w(0,y) > 0. At negative infinity we 
had exponential convergence \pi(x,y) — 5\ < Ce^ x , and we chose the supersolution 
z = e Xx to decay slowly (A > was chosen small enough, for example A = Co/2m8), 
hence \z\ = \pi — p 2 1 < Ce^ x <C \z\ and w(—oo,y) = 0. The Minimum Principle 
applied to L[w] = finally shows that w(x, y) > for x < 0, and therefore p\ > p 2 on 
the whole cylinder D = K x T 1 if p 2 is slided far enough to the right. 

Slowly sliding back to the left we obtain a first critical translation p 2 , after which we 
cannot keep translating to the left without breaking p\ > p 2 (this critical translation 
exists because sliding p 2 far enough to the left the two solutions must cross at some 
point). By continuity we have that z* = p\ — p* 2 > 0, and we claim that there exists 
a contact point (x ,y ) 6 D such that z*(x ,y ) = 0. Temporarily admitting this, we 



obtain a contradiction as follows: z* > satisfies (6.1 ), which is uniformly elliptic with 
bounded zero-th order coefficient, and attains a minimum point z* = in D. The 
Minimum Principle shows that z* = 0, thus contradicting q^i > q 2y i z(+oo,y) = 
+oo. 

In order to obtain such a contact point, assume by contradiction that p\ > p* 2 on 
the whole cylinder: condition g^i > g 2 ,i shows that p\ — p\ > C a > on any sub- 
cylinder x > —a for any a > large and some constant C a . We may then slide p 2 
slightly further to the left in such a way that pi > p 2 if x > a. Repeating the above 
comparison argument for x < a, we see that p\ > p 2 also holds to the left, hence on 
the whole cylinder. This contradicts the fact that p 2 was a critical translation. 

We just proved that qi i > g 2 ,i cannot hold, and by symmetry pi -h- p 2 we obtain 

We may now repeat the very same argument to show that q 12 = q 2t2 , and so forth {q^k — 
q2,k =^ Qi,k+l — Q2,k+i)- n 

We can now prove uniqueness of the 5-solutions: 



Proof, (of Theorem 6.1) Let p\,p 2 be two 5-solutions; we pin as before p\ once and for all, 
and only slide p 2 . Let us stress that both solutions have now the same coefficients q\...q^ in 
their asymptotic expansion at infinity 

i = l,2 x — > +oo : Pi(x, y) = cx + x (^qix^™ + ... + q N x~^j + q* + 0^(1), 
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except maybe for the last two terms (the lower order q* + Oj(l)). We recall that z := p\ — p 2 



satisfies (6.1 ) of the form L[z] = 0, and remark the following: for any r-translation p 2 (x—r, y), 



uniqueness of the coefficients q\...qN shows that 

Pi(x,y) - p 2 (x - r,y) = ct + g* - q 2 + o(l). (6.2) 

+ 0O 

This means that, depending on the translation, only two scenarios are possible at infinity: 
either [pi — p 2 )(+oo,y) = 0, either [p x — p 2 ](+oo,y) = est ^ 0. 

We showed previously that sliding p 2 far enough to the right p± > p 2 must hold, and that 
slowly sliding back to the left there exists a first critical translation p 2 such that p\ >p 2 - 
Similarly translating p 2 far enough to the left we have that p\ < p 2 , and there exists a first 
critical translation p 2 coming from the left side such that p\ < p 2 . 

1. If there exists a contact point pi(xo,yo) = P 2 (%o,yo) then z = p\ — p 2 is nonnegative 
(because p 2 is a critical translation coming from the right side), satisfies an elliptic 
equation C[z] = with bounded zero-th order coefficient, and attains an interior 
minimum point z(xo, yo) = 0. The classical Minimum Principle shows that z = 0, 
meaning precisely that p\ can be deduced from p 2 by translation. 
We may therefore assume that no such contact point exists, and the only possible 
scenario is therefore that [pi — p 2 ](+oo,y) = (otherwise [p\ — p 2 ](+oo,y) = est > 



according to (6.2) and we could slide p 2 a little further to the left as in the proof of 



Proposition 6.2, thus contradicting the fact that p 2 is critical). 



2. Similarly arguing for z, we may assume that [pi — p 2 ](+oo,y) = 0. 

3. As a consequence \p 2 — p 2 ](+oo,y) = 0, and therefore p 2 = P 2 - We conclude recalling 
that we constructed p 2 < pi < p 2 , hence p\ =p 2 = P 2 - 

□ 
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